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Abstract. The dynamical correlation functions in the one-dimensional electronic 
systems show power-law behavior at low energies and momenta close to integer 
multiples of the charge and spin Fermi momenta. These systems are usually referred 
to as Tomonaga-Luttinger liquids. However, near well-dfined lines of the {k,Lu) plane 
the power-law behaviour extends beyond the low-energy cases mentioned above, and 
also appears at higher energies leading to singular features in the photoemission 
spectra and other dynamical correlation functions. The general spectral-function 
expressions derived in this paper were used in recent theoretical studies of the finite- 
energy singular features in photoemission of the organic compound tetrathiafulvalene- 
tetracyanoquinodimethane (TTF-TCNQ) metallic phase. They are based on a so called 
pseudofermion dynamical theory (PDT), which allows us to systematically enumerate 
and describe the excitations in the Hubbard model starting from the Bethe-Ansatz, as 
well as to calculate the charge and spin objects phase shifts appearing as exponents 
of the power laws. In particular, we concentrate on the spin-density m — > limit 
and in effects of the vicinity of the singular border lines, as well as close to half 
filling. Our studies take into account spectral contributions from types of microscopic 
processes that do not occur for finite values of the spin density. In addition, the specific 
processes involved in the spectral features of TTF-TCNQ are studied. Our results are 
useful for the further understanding of the unusual spectral properties observed in 
low-dimensional organic metals and also provide expressions for the one- and two- 
atom spectral functions of a correlated quantum system of ultracold fermionic atoms 
in a ID optical lattice with on-site two-atom repulsion. 
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1. INTRODUCTION 

The low-energy physics of correlated one- dimensional (ID) problems has some universal 
properties described by the Tomonaga-Luttinger liquid (TLL) [Jj. In turn, a 
pseudofermion dynamical theory (PDT) beyond the TLL [2], |3] was recently used to 
study the finite-energy singular features in photoemission of the organic compound 
tetrathiafulvalene-tetracyanoquinodimethane (TTF-TCNQ) metallic phase While 
the PDT was originally introduced for the ID Hubbard model, more recently other 
methods for the study of finite-energy spectral and dynamical functions of ID correlated 
systems were introduced [5], E]. Both the finite-energy spectral- weight distributions 
studied by the PDT for the ID Hubbard model and the methods of Refs. [5|.6J for other 
ID correlated problems include power-law singularities near well-defined branch lines 
with exponents depending on the interaction strength and the excitation momentum. 
In the limit of low energy the PDT correlation- and spectral-function expressions recover 
the usual low-energy TLL results [7j 

Besides a renewed interest in the unusual spectral and dynamical properties of 
quasi-lD organic compounds [8], [9], recently there has been an increasing interest in 
those of interacting ultracold fermionic atoms in ID optical lattices [TO]. Quantum 
effects are strongest at low dimensionality leading to unusual phenomena such as 
charge-spin separation at all energies P [9]. Thus, the further understanding of 
the microscopic mechanisms behind the unusual spectral properties observed in low- 
dimensional correlated systems and materials is a topic of high scientific interest. 

The ID Hubbard model is one of the few realistic models for correlated electrons 
in a discrete lattice for which one can exactly calculate all the energy eigenstates and 
their energies pTl [T2| [13] . It includes a first-neighbor transfer-integral t, for electron 
hopping along the chain, and an effective on-site Coulomb repulsion U. For finite-energy 
values the metallic phase of this model goes beyond the low energy behavior described 
by the usual TLL [1] and thus the study of spectral functions is a very involved many- 
electron problem. Fortunately, the construction of a pseudofermion description by means 
of a unitary transformation which slightly shifts the discrete momentum values of the 
corresponding pseudoparticles of Refs. [TUITS] leads to the PDT, whose energy spectrum 
has no residual-interaction energy terms. Therefore, such a description is suitable for 
the derivation of explicit expressions for these functions [21 [3l [7]. 

The PDT is a generalization for finite values of U of the scheme introduced in Refs. 
[T6| [T7] for large U values. It profits from the use of a description of the exact energy 
eigenstates in terms of occupancy configurations of several branches of pseudofermions. 
The ground state has finite occupancy of charge c and spin si pseudofermions only. 
Under the ground-state - excited-energy-eigenstate transitions, the pseudofermions and 
pseudofermion holes undergo elementary scattering events with the pseudofermions and 
pseudofermion holes created in these transitions. This leads to excited-state, interaction, 
density, and two-momenta dependent two-pseudofermion phase shifts. The point is that 
the one- and two-electron spectral functions can be expressed in terms of pseudofermion 
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determinants which are a functional of such phase shifts. Use of the PDT reveals that 
for finite values of U/t all singular one-electron spectral features of the model are of 
power-low type, controlled by negative exponents. Furthermore, the PDT line shapes 
associated with such exponents were found to correspond to the unusual charge and 
spin spectral features observed by photoemission experiments for the whole finite-energy 
band width in quasi-lD organic metals [18]. (The use of the Dynamical Density Matrix 
Renormalization Group method [19] leads to results consistent with those obtained by 
the PDT.) Furthermore, when combined with the Renormalization Group, the use of 
the PDT reveals that a system of weakly coupled Hubbard chains is suitable for the 
successful description of the phase diagram observed in quasi-lD doped Mott-Hubbard 
insulators [20] . 

The low-energy physics of the model corresponds to the universal TLL behaviour 
and was studied by different techniques, such as conformal-field theory [21] and 
bosonization [22]. There are many investigations where the low-energy conformal 
invariance was combined with the model exact Bethe-ansatz solution in the study of 
the asymptotics of correlation functions and related quantities [H, [22] . As mentioned 
above, the studies of Ref. [7] confirm that in the limit of low energy the general finite- 
energy spectral and correlation function expressions provided by the PDT in Ref. [2] 
recover the correct behaviour given by conformal-field theory. 

The main goal of this paper is to provide the details of an extension of the PDT 
introduced in Refs. [2], [3l [7] to initial ground states with spin density m ^ which was 
used in recent theoretical studies of the finite-energy singular features in photoemission 
of the TTF-TCNQ metallic phase [1]. The point is that such an extension involves 
spectral contributions that do not occur for finite values of the spin density and hence 
are not taken into account by the expressions of Refs. [21 [3[ [7] . In order to introduce the 
generalized spectral-function expressions used in the recent studies of the TTF-TCNQ 
spectral features presented in short form in Ref. [1] we present here such an extension 
of the PDT. 

The unusual spectral properties are mainly determined by the occupancy 
configurations in the excited states of the two pseudofermion branches that have finite 
occupancies for the ground state. Those are the charge cO pseudofermion and spin si 
pseudofermion branches. The ground-state corresponds to charge cO pseudofermion (and 
spin si pseudofermion) finite occupancy for canonical- momentum values in the range 
1^1 < Q'fco — (and |g| < qp^^ = kpi) and charge cO pseudofermion-hole (and spin 
si pseudofermion-hole) finite occupancy for canonical- momentum values in the range 
2kp < \q\ < = vr (and kpi < \q\ < q^si = kp])- The values of the four deviations 
Agi^co, ±1 and /S.qpsi^±i in the two charge ±5'^co ~ '^'2,kp and two spin ig^^i = i^Fj, 
Fermi points under each ground-state - excited-energy-eigenstate transition play a major 
role in the PDT. 

We denote such deviations by Aqpau,L where au = cO, si and l = ±1. The 
spectral-weight distributions are controlled by the values of the following associated 
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four parameters, 

2A^.= ; au = cO,sl; . = ±1, (1) 

where L >> 1 is the ID lattice length within the use of periodic boundary conditions. 
According to the PDT, the contributions from ground-state transitions to subspaces 
spanned by sets of excited energy eigenstates with the same values for the two charge 
parameters 2A^^ and two spin parameters 2A^j^^ fully determine the momentum k and 
energy u dependence of the general finite-energy spectral functions in the small {k, uj)- 
plane region associated with the energy and momentum spectrum of these excited states 

EE]. 

The general PDT finite-energy spectral-function expressions given in Ref. [2] refer 
to the metallic phase for initial ground states with spin densities m in the range 
< m < n, where n is the electronic density such that < n < 1. Actually, the 
expressions found in Ref. [2] refer to values of n such that the Fermz-point charge 
velocity Vco is larger than the Fermi-point spin velocity Vgi- For finite values of U/t, 
this excludes densities n in the vicinity of half filling. One of our goals is to extend the 
studies of that reference to electronic densities such that Vsi > fco- 

The general expressions of the exponents that control the singular features of the 
above finite-energy spectral function expressions provide the correct zero spin-density 
values in the limit m — 0. Moreover, as given in Eq. (55) of Ref. [2], for excitations 
such that 2A^^ ^ for the two au = cO charge and two au = si spin parameters, 
the corresponding four pseudofermion relative weights have the asymptotic expression 
provided in that equation. It follows that for zero-spin excitations such that the four 
parameters 2A^^^ where au = cO, si and l = ±1 are finite, the general convolution 
function and corresponding pre-factor function Fo{z) given in Eqs. (61) and (62) of Ref. 
[2], respectively, provide the correct zero spin-density contributions to the spectral- 
function expressions in the limit of zero spin density. The point is that for densities 
in the ranges < n < 1 and < m < n all ground-state - excited-energy-eigenstate 
transitions lead to finite values for these four parameters. 

While some m ^ one- and two-electron excitations also lead to finite values for 
the two au = cO charge and two au = si spin parameters 2AJ;^^, there are also m — 
excitations for which one (or both) spin parameter (s) 2Af^ vanishes (or vanish). In 
that case the corresponding si, t pseudofermion relative weights have not the asymptotic 
expression provided in Eq. (55) of Ref. [2]. Nonetheless, we find in this paper that the 
contributions to the spectral-function expressions from ground-state transitions to the 
excited energy eigenstates which span such excitations lead to convolution functions of 
the same general form as that provided in Eq. (61) of Ref. [2]. The only difference is 
that the pre-factor Fq{z) of such convolution functions given in Eq. (62) of that reference 
is replaced by another suitable function derived in this paper and given in Appendix A. 

Since these convolution functions fully control all the different spectral-function 
contributions given in Eqs. (66), (68), and (70) of Ref. [2], the derivation of all m = 
contributions requires the use of such suitable pre-factors, which were taken into account 
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in the studies of the TTF-TCNQ spectral features of Ref. [1] and we calculate in this 
paper. Interestingly, the different expressions found here for the pre-factor Fq{z) are 
such that the corresponding pre-factors of the space and time asymptotic expressions of 
correlation functions |7j| are continuous functions of m as m — >■ 0. 

The extension of the expressions for the one- and two-electron spectral-weight 
distributions introduced in Ref. [2] for the ID Hubbard model to the regime where 
m —>■ involves taking into account contributions from transitions to subspaces spanned 
by excited states such that one, two, three, or even four out of the four parameters 2A^^ 
where au = cO, si and l = ±1 vanish. Note that according to Eq. ([1]), the value of the 
au, L pseudofermion canonical-momentum Fermi points associated with 2A^^ = values 
remains unchanged under the corresponding ground-state - excited-state transitions. 
Such contributions do not exist for the m > regime addressed in Ref. [2] but must 
be taken into account in the quantitative study of the spectral-weight distributions of 
the correlated metal at zero spin density. In this paper we also extend the one- and 
two-electron spectral- weight distributions of Ref. p] to all electronic densities of the 
metallic phase. Although we do not consider the half-filling Mott-Hubbard insulator 
such that Vco = 0, our expressions refer to all electronic densities of the metallic phase 
and thus also for those in the vicinity of 1. 

Another problem solved in this paper is the derivation of explicit expressions for the 
spectral functions in the vicinity of the singular border lines, which again were used in 
the recent studies of the TTF-TCNQ spectral features of Ref. [1]. For that one-electron 
problem such border lines correspond to processes for which the extra charge and spin 
objects created upon removal or addition of the electron have exactly the same velocity. 
Such an equahty of the charge and spin velocities occurs at well-defined lines in the (fc, u) 
plane and leads to power-law singular behavior along such a border lines. In contrast to 
the power-law branch-line singular features studied in Ref. [2], which are controlled by 
momentum, U/t, and density dependent negative exponents, the border-line power-law 
singular features are controlled by an universal exponent given by —1/2, as given in 
Eq. (2) of Ref. [1]. Further details on the processes involved in the applications to 
TTF-TCNQ are also reported. 

The paper is organized as follows: In Sec. [2] we introduce the model and provide 
basic information about the pseudofermion description needed for our studies. In Sec. 
[3] the general expressions required for the study of the finite-energy spectral-weight 
distributions of the metallic phase for initial m ground states are calculated. That 
includes derivation of explicit expressions of the spectral functions in the vicinity of the 
singular border lines for density ranges < n < 1 and < m < n. In Sec. |4]we use the 
expressions obtained in the previous section and in Ref. |2j to provide further details on 
the processes that contribute to the unusual photoemission spectrum of TTF-TCNQ. 
Finally, the concluding remarks are presented in Sec. [5l 
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2. THE PROBLEM AND THE PSEUDOFERMION DESCRIPTION 

Our study focus on finite-cj A/'-electron spectral-weight distributions of the following 
general form, 

B'^ik,u;) = J2\{f\0'Afik)\GS)\'5{uj-l[Ej-EGs]); Ito > , (2) 
/ 

where / = ±1. We focus our attention on the cases of more physical interest which 
correspond to A/" = 1, 2. In the above expression the general A/'-electron operators 
0^{k) = 0\j-{k) and Ojl{k) = Oj\f{k) carry momentum k, the / summation runs over 
the excited energy eigenstates, the energy Ef corresponds to these states, and Eqs is 
the initial ground-state energy. The local operator 0^,j = 0'^j\r,3 "-"^ ^JV'lj = ^■N',] 
related to the corresponding momentum-representation operator 0^j^{k) of Eq. ([2]) by 
a Fourier transform. As in Ref. [2], we use in expression ([2]) a momentum extended 
scheme such that k G {—oo, +oo). 

We consider weight distributions ([2]) that refer to the Hubbard model in a ID 
lattice with periodic boundary conditions and units such that the Planck constant and 
electronic lattice constant are one, 

H = -t^[c]^^Cj+i,^ + h.c] + UY^ %T^i,i • (3) 

j, o- j 

Here cj ^ (and Cj^^) creates (and annihilates) one spin-projection a =], | electron at site 

' t 

j = 1, 2, Aq and fij^cr = c]-,^ Cj^a- Let A = A^ + A^ be the electronic number, Aq the 
number of lattice sites, and n„ = N„/L = N„/Na- is assumed to be even and very 
large. The electronic densities n = + and spin densities m = — rij are in the 
ranges < n < 1 and < m < n, respectively. Except for corrections of order of l/L, 
the Fermi momenta are given by kp = 7m/2 and kp^ = 7m„. 

The concept of a rotated electron plays a key role in the pseudofermion description. 
Concerning its relation to the holons, spinous, and cO pseudoparticles whose occupancy 
configurations describe the energy eigenstates of the model ([3]), see Ref. [Hj. Our 
studies do not involve directly the holons and spinous as defined in that reference. 
The charge cu pseudofermions (and spin su pseudofermions) such that z/ = 1,2,3,... 
are 2z/-holon (and 2z/-spinon) composite quantum objects whose discrete momentum 
values are slightly shifted relative to those of the corresponding cu pseudoparticles 
(and spin su pseudoparticles) studied in Refs. [Tll[l5]. Such momentum shifts cancel 
exactly the residual-interaction terms of the pseudoparticle energy spectrum. Otherwise, 
pseudoparticles and pseudofermions have the same properties. 

According to the PDT of Refs. [2], [7], the charge cO pseudofermions and spin si 
pseudofermions play the major role in the spectral properties. The holons (and spinous) 
which are not part of 2z/-holon composite ciy pseudofermions (and 2z/-spinon composite 
su pseudofermions) are the Yang holons (and HL spinous). Those are invariant under 
the electron - rotated-electron unitary transformation and hence have a non-interacting 
character and do not contribute to the matrix elements between energy eigenstates 
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of the spectral-weight distributions studied in this paper. We denote the numbers of 
au pseudofermions and au pseudofermion holes by Nau and N^^^, respectively, where 
a = c, s and z/ = 0, 1, 2, ... for a = c and z/ = 1, 2, ... for a = s. (The value of N^^ is 
given in Eqs. (B7) and (B8) of Ref. [H].) 

Alike in Ref. [2], we use in this paper the notation au ^ cO, si branches, which 
refers to all av branches except the cO and si branches. Moreover, the summations 
(and products) ^ 

ai^=cO, si' and X!]ai/:^cO, si 

(and n„., n 1' a^*^ riai/^cCsl) 
run over all au branches with finite au pseudofermion occupancy in the corresponding 

state or subspace, the cO and si branches only, and all au branches with finite au 

pseudofermion occupancy in the corresponding state or subspace except the cO and si 

branches, respectively. 

The pseudofermion description refers to a Hilbert subspace called the pseudofermion 

subspace (PS) in Ref. in which the jV-electron excitations 0\j-{k)\GS) are contained. 

The PS is spanned by the initial ground state and the excited energy eigenstates 

originated from it by creation, annihilation, and particle-hole processes involving the 

generation of a finite number of active pseudofermion scattering centers, Yang holons, 

and HL spinous plus a vanishing or small density of low-energy and small-momentum 

au = cO, si pseudofermion particle-hole processes. It is convenient to classify these 

processes into three types, called processes (A), (B), and (C), as further discussed in 

the ensuing section. 

The a //-pseudofermion discrete canonical-momentum values have a functional 
character and read, qj = qj + Q^^{qj)/L = [2n / L]!""" + Q^^{qj)/L where j = 1, 2, iV*^, 
and N*j^ = N^u + N'^y Here is a ^(gj)/2 is a av pseudofermion scattering phase 
shift given by, 

N", , 

a' u' 

QL{Qj)/2 = Y1 E ^'^au,a'u'{q,,qr)ANo,,,,{q,,)- j = 1,2, , (4) 

a'u' j'=l 

where ANau{qj) = AAfau{qj) is the bare-momentum distribution function deviation 
A.Nav{qj) = Nau{qj) — N^uilj) corresponding to the excited energy eigenstate. This 
deviation is expressed in terms of the bare-momentum qj = [2itI°"^]/ L, which is carried 
by the au pseudoparticles, where 1°"^ are the quantum numbers provided by the Bethe- 
ansatz solution jl4j . 

Although the au pseudoparticles carry bare- momentum qj, one can also label the 
corresponding au pseudofermions by such a bare-momentum. When we refer to the 
pseudofermion bare- momentum qj, we mean that qj is the bare-momentum value that 
corresponds to the pseudofermion canonical momentum qj = qj + Q^j^{qj) / L. For the 
ground state the pseudofermion numbers are given by A^^co = A^, A^si = N^, N^y = for 
au 7^ cO, si. We call A'^q and N^-^ the ground-state cO and si pseudofermion numbers, 
respectively. As mentioned in the previous section, the ground-state au = cO, si bare- 
momentum distribution functions are such that there is pseudofermion occupancy for 
IqI < Qpau unoccupancy for q^^^ < \q\ < qau^ where in the thermodynamic hmit 
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the Fermi-point values are given by, 

Ifco = 2fcF ; qpsi = kpi ■ (5) 

Moreover, for that state the hmiting bare-momentum values of both the av = cO, si 
and au ^ cO, si bands read, 

9co = ^ i = ; (fcv = [7r-2A;j7] , > ; = [kpi-kFi] , z/ > 1 .(6) 

The ground-state au = cO, si densely-packed bare-momentum distribution functions 
N^^iq^) are given in Eqs. (C.1)-(C.3) of Ref. [g. 

Under the ground-state - excited-energy-eigenstate transitions, the au pseud- 
ofermions and au pseudofermion holes undergo elementary scattering events with the 
a'u' pseudofermions and a'u' pseudofermion holes created in these transitions [2]. This 
leads to the elementary two-pseudofermion phase shifts n ^au,a'u'{,(lj, Qj') on the right- 
hand side of the overall scattering phase shift (jl]), as further discussed in Ref. |23] . 
Moreover, within the PDT the overall au pseudofermion or hole phase shift, 

Qa.fe)/2 = gL/2 + QLfe)/2, (7) 

controls the spectral weight distributions. Here Qau{<lj)/L gives the shift in the discrete 
canonical-momentum value qj that arises due to the transition from the ground state to 
an excited energy eigenstate and Q^uf^ can have the values Qa^f^ = 0, ±tt/2 |2l [23] . 
In this paper we use boundary conditions such that Q^^yjl = 0, —sgn{k) tt/2, where 
k is the excited-state momentum relative to that of the initial ground state. Here we 
assume that for the latter state N/2 and N are odd and even numbers, respectively. 
Q^^,/L gives the shift in the discrete bare-momentum value qj that arises as a result of 
the same transition. 

The au pseudofermion creation and annihilation operators f^.^ai^ and fq^^au, 
respectively, have exotic anticommutation relations [2]. These anticommutators involve 
the overall phase shifts ([7]) and play a key role in the spectral properties. There are 
corresponding local operators fl^^au fxj,au- Here Xj where j = 1,2, ...,N*j^ are the 
spatial coordinates of a au effective lattice with N*^ sites. Each of the Nau occupied 
sites of such a lattice correspond to well-defined occupancy configurations of 2i/ sites 
of the rotated-electron lattice. It turns out that the operators f^^^au fxj,au have 
simple anticommutation relations. Indeed, the exotic anticommutation relations of the 
operators and f^^^au and fq^^au result from the form of the discrete canonical momentum 
values q, = g, +QM/L = [27r/L]/- + QM/L. 

There is a simple relation between the f^^^co a'^'^ fxj,co operators and those of 
the rotated electrons such that the former operators have anticommutation relations. 
To show that for branches au ^ cO the operators j\^^av ^^cl fx^,av also have 
anticommutation relations is a much more involved problem. The composite av ^ cO 
pseudofermions emerge from from corresponding av ^ cO bond particles through a 
suitable Jordan- Wigner transformation [21]. The point is that such av ^ cO bond 
particles live on the corresponding av effective lattice and behave there as hard-core 
bosons. That interesting problem will be studied in detail elsewhere. 
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The charge and spin pseudofermion Fermz-point group velocities Vco and Vgi also 
play an important role in our study. The velocity Vau with au = cO, si is a particular 
case of the momentum- dependent group velocity Vau{<i)- Such velocities appear in all 
the spectral-weight distribution expressions of the metallic phase and are given by, 

Vauiq) = — TT — , all branches ; v^u = Vauiqpau) > = cO, si . (8) 
oq 

In the first expression eauiq) stands for the au pseudofermion energy dispersion defined 
by Eqs. (18)-(20) of Ref. [2j. In this paper and its Appendix A we use often a convention 
according to which the az/ = cO, si branch is that whose pseudofermion Fermz-point 
group velocity Vap is such that, 

Vau = T^^T^{VcO, Vsl} . (9) 



3. SPECTRAL- WEIGHT DISTRIBUTIONS FOR THE METALLIC 
PHASE AND 0<m<n 

3.1. General spectral-function expressions 

The pseudofermion elementary processes that generate the PS from the initial ground 
state belong to three types: 

(A) Finite-energy and finite-momentum elementary cO and si pseudofermion 
processes away from the corresponding Fermi points involving creation or annihilation 
of a finite number of pseudofermions plus creation of av 7^ cO, si pseudofermions with 
bare-momentum values different from the limiting bare- momentum values ig'^j,; 

(B) Zero-energy and finite-momentum processes that change the number of cO 
and si pseudofermions at the l = sgn(g) 1 = +1 right and l = sgn(g) 1 = — 1 left 
cO and si Fermi points - these processes transform the ground-state densely packed 
bare-momentum occupancy configuration into an excited-state densely packed bare- 
momentum occupancy configuration. Furthermore, creation of a finite number of 
independent —1/2 holons and independent —1/2 spinous, including —1/2 Yang ho Ions, 
—1/2 HL spinous, and —1/2 holons and —1/2 spinous associated with cu pseudofermions 
of limiting bare momentum q = ±g^^, = ±[7r — 2kF] and su pseudofermions of limiting 
bare momentum q = ±g°^, = ±[kp-^ — kp^], respectively; 

(C) Low-energy and small-momentum elementary cO and si pseudofermion particle- 
hole processes in the vicinity of the l = sgn(q') 1 = +1 right and l = sgn(g) 1 = —1 left 
cO and si Fermi points, relative to the excited-state au = cO, si pseudofermion densely 
packed bare-momentum occupancy configurations generated by the above elementary 
processes (B). 

Such processes generate excitations which can be classified by the values of a 
set of numbers and number deviations. For instance, N^'^^ is the number of finite- 
momentum and finite-energy av = cO, si pseudofermion particle-hole processes of type 
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(A). The quantum number l = sgn(g)l = ±1 refers to the right pseudfermion movers 
{l = +1) and left pseudfermion movers {l = —1) and AN^^^ ^ such that AN^^ is the 
deviation in the number of au pseudofermions at the right (+1) and left (—1) Fermi 
points generated by the elementary processes (B). In turn, the deviation in the number 
of au = cO, si pseudofermions created or annihilated away from these points by the 
elementary processes (A) is denoted by AN^^f . 

The actual number of au pseudofermions created or annihilated at the right (+1) 
and left (—1) Fermi points by the processes (B) is denoted by AN^l'f.^^. It is such 
that AN^j^ ^ = AN^'jf^ + lQ^^^/2ti, where Q^^vl'^ scattering-less phase shift on the 

right-hand side of Eq. ([71). Furthermore, ^ refers to the av ^ cO, si branches and 
is the number of av pseudofermions of limiting bare momentum q = cq^^ such that 
L = ±1 created by the elementary processes (B). The number of au pseudofermions 
created away from the limiting bare-momentum values by the processes (A) is called 

av ■ 

A ground-state - excited-energy-eigenstate transition generated by elementary 
processes (A) and (B) leads to the energy and momentum spectrum lAE and ZAP 
given in Eqs (28) and (29) of Ref. [2], respectively, where / = +1 or / = — 1 depending 
on the specific spectral function (j2]) under consideration. Each of such excited states is 
associated with a well defined point {lAE ,IAP) in the (/c, tu) -plane. A key property of 
the PDT is that the set or tower of excited states generated by the elementary processes 
(C) from each excited energy eigenstate generated by the processes (A) and (B) have 
the same values for the two charge parameters 2A^^ and two spin parameters 2A^j^^ of 
Eq. ([1]) as the latter state. The transitions to that tower of excited states generates 
the spectral weight in the vicinity of the corresponding point {lAE^lAP). Thus, that 
weight is associated with the same value of the following functional, which plays a major 
role in the spectral properties. 



where av = cO, si and t = ±1. Note that Co equals the sum of the four parameters of 
Eq. (lH), which can be expressed in terms of the overall scattering phase shift of Eq. (jl]), 
as given in the second expression of Eq. (ITOl) . Such parameters are functionals of the 
pseudofermion occupancy configurations which describe the excited states generated by 
the elementary processes (A) and (B). 

The corresponding functional expressions are given in Eqs. (12) and (40) of Ref. 
[2]. Thus, for each excited state generated from the initial ground state by the processes 
(A) and (B) there is a subspace spanned by the set of excited states generated by the 
processes (C) from the former excited state. Given the linear av = cO, si pseudofermion 
energy dispersion near the Fermz points, the processes (C) lead to small momentum and 



Co 



2A+^ + 2A,-oi + 2A+i + 2A 



-1 

■si 5 




(10) 
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energy values such that, 

k' = ^ L ^ m^^, I, ; uj' = ^ '^Vaumo,u,c- (H) 

01^=00, si t=±l «!^ = cO, si t=±l 

Here mQ.^, ^ is the number of elementary au pseudofermion particle-hole processes of 
momentum l[27t/L]. Thus, the elementary processes (C) generate a set of excited energy 
eigenstates with energy lAE and momentum ZAP given by those of the initial excited 
state state generated from the ground state by the processes (A) and (B), plus the small 
energy u' and momentum k' provided in Eq. ( ITTi) . 

Within the PDT the general A/"-electron spectral functions ([2]) factorize in terms of 
au pseudofermion spectral functions. The probability amplitude Aaif^ appears in the 
expressions of the latter functions for au = cO, si. It is associated with the canonical- 
momentum densely packed configurations generated by the elementary processes (A) 
and (B) after the matrix elements between the ground state and the excited states 
generated by these processes are computed. Such a probability amplitude corresponds 
to the au = cO, si pseudofermion spectral-function lowest-peak weight given in Eq. (48) 
of Ref. It has the following approximate behaviour. 



4°;°N7 t'^W^' «^ = cO,s1; . = ±1. (12) 



Here 2A^^ is the functional given in Eq. ([T]), and fau,i. reads. 



fau,, = \j f(Qau{iq%av) + sgn(A;)7rj ; = JJ ; au = cO, si . (13) 

In this expression, k stands for the excited-state momentum relative to the initial ground 
state, f{Q) = f{—Q) is the function defined in Ref. [IT], which occurs on the right- 
hand side of Eq. (24) of that reference, and fau appears in spectral-function expressions 
introduced below. Moreover, 5'°^, is a n, m, and U/t dependent constant such that 
^co S^si ~^ 1 both for U/t ^ and for U/t ^ oo and m — » 0. (From Ref. [TT] we learn 
that S^Q sin(7m) for U/t ^ oo and m 0, and thus S^^ ^1/ sin(7m) in such a 
limit.) It is useful to introduce the following quantity, 

D = TT '^Q^ f (i4\ 

11 (CO M2A+i+2A-,i] ■ ^ ' 

au=cO,sl y'^au) 

Another important peace of the au = cO, si pseudofermion spectral functions is the 
relative weight aau,L{^au,L) generated by the elementary processes (C) in the vicinity 
of the right {l = +1) and left [t = —1) Fermi points. When 2A''^^ > the weight 
o-av,i{jnav,i) and its asymptotic expression are given in Eqs. (52) and (55) of Ref. p], 
respectively. The point is that when the excited states generated by the processes (A) 
and (B) are such that = the relative weight aau,L{^au,L) reads instead, 

aau,,ifnau,,) = Sm^^^^o; 2A^^ = 0; au = cO,sl; l = ±1 . (15) 
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Since for densities in the range < n < 1 and < m < n the two charge parameters 
2A^^ and two spin parameters 2A^^^ are finite for all excited states, this case was not 
considered in the studies of Ref . [2] . 

As further discussed below, the general A/'-electron spectral function given in Eq. 
([2]) can be expressed in terms of the charge cO and spin si spectral functions provided 
in Eq. (56) of Ref. [2]. Here we express the latter functions as follows, 

B'lAk', u') = ^ dk" du" B^:^ {k' - k\ J - J' 

X B'£'f (^k",uj"y, au = cO,sl; l = ±1 . (16) 

The studies of Ref. [2] refer to the case where the two charge parameters 2A^^ and two 
spin parameters 2Af^ are finite. In that case the four relative weights aau, i.i'mau, l) such 
that au = cO, si and l = ±1 and their asymptotic expressions are of the form provided 
in Eqs. (52) and (55) of Ref. [2], respectively, and thus the functions BQ^^{k' , u') on 
the right-hand side of Eq. ( fT6l) are given by [2] , 

4(0,0) , / / 

BQlSk.u) = "-■4 2vrt; In ) ~ ~^ 

'-'au ^ ^" '-'qj^/ ' a ' '^au 

~ ^--^ ( Y^^--'5(k' - (17) 

where av = cO, si and l = ±1. The second expression corresponds to the asymptotic 
behaviour valid for small finite values of lu'. 

However, when 2A^^ = the corresponding relative weight aau, iii^au, l) is of the 
form provided in Eq. ( |T5l) and thus the spectral function BQ^^{k',u') instead of being 
given by Eq. (ITTj) reads, 

B^-:lik',uj') = ^A^^;':i6ik')6i^') 

/ 90 

-'2nU,^J^6{k')6{uj'); au = cO,sl; t = ±1 . (18) 

In order to reach the correct spectral-function expressions for all electronic densities 
of the metallic phase and for m — > 0, several cases must be considered. The first 
corresponds to excited states such that the two charge parameters 2A^^ and two spin 
parameters 2Af^ are finite. The two functions on the right-hand side of Eq. (fT6|) are of 
the form given in Eq. (fT7|) and thus we find, 

^ ^(00) TT flluj' + LVauk']- 



4:7CVau/Na 



fau -TT 9(/[^^ + LVa„ k']) d[uj' + LVg^ fc^] \ 2A;..-1 



4vrt;,,^o, Ai r(2A^J KAttv^^S^, 



au '-'au 



where au = cO, si. The second expression corresponds to the asymptotic behaviour 
valid for small finite values of l[u!' + LVau k']. This function is that already given in Eq. 
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Figure 1. Full theoretical line shape of the one-electron removal spectral- weight 
distribution found in Ref. i4J to fit the corresponding spectral features of TTF-TCNQ. 
The spectrum shown here is rotated relative to that shown in Fig. 2 of that reference. 
It includes both the TTF related spectral features for n — 1.41; t = 0.35 eV; U/t = 5.61 
and those of TCNQ for n = 0.59; t = 0.40 eV; U/t = 4.90, respectively. 



(56) of Ref. |2]. (See the note [25] concerning a discrepancy of a factor 2 between the 
function given here and that provided in Eq. (56) of Ref. [2] and how a misprint in the 
latter equation becomes an error, by propagating to other expressions of Refs. [21 [7].) 

The second case occurs when the excited states are such that for the functions on 
the right-hand side of Eq. f|T6|) one has that 2A^j^ > and 2A~^ = 0. Then one must 
use in that equation the expression (JTTj) for the function Bg^^^k' ,u') and that given in 
Eq. ( ITSl) for the function BQ~''f{k',uj'). This leads to the following expression for the 
function on the left-hand side of Eq. (JT6l) . 



fau 



A, 



(0,0) 



luj' 



where av = cO. 



^«.r(2AL, 

si and L = ±1. 



Q{liu') 



lu' 



2lT Vau S^u 



(20) 



Finally, the third case refers to excited states such that for the functions on the 
right-hand side of Eq. f|T6|) one has that 2A^j^ = 2A~^ = 0. In this case one must use in 
that equation the expression ( fTSl) for both the functions BQ^^{k',uj') and BQ~''j^{k',u!'). 
One then reaches the following expression for the function on the left-hand side of Eq. 



N, 
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Figure 2. The region of the {k,uj) plane with a finite one-electron removal spectral 
weight obtained by running over all the momentum values of the cO and si momentum 
distribution deviations of Eqs. ([5^ and (|33p for U/t ^ 100, n = 0.59, and m 0. 
The corresponding branch and border lines are also plotted. 




kp Skp n Q Skjp n 



Figure 3. The same region of the (fc, lS) plane and branch and border lines as in Fig. 
2 for Ujt = 4.9, n = 0.59, and m ^ 0. 



^lT:f^,SlJ{U)b{J)- az/ = cO, Si. (21) 

The general PDT expression for the A/'-electron spectral function given in Eq. ([2]) 
reads [2], 

oo 

where Cg = 1, / = ±1, and for densities < n < 1 the summation over the numbers 
N^^^ is limited to finite values. On the right-hand side of Eq. (!22|) . c\ is the constant 
of the operator expressions given in Eqs. (32)- (34) of Ref. (3j such that c- — > as 
f//t— >ooforz>0 and the function i?''*(/c, uj) is defined in Eq. (44) of Ref. [2J. The 
latter function is fully defined by the related function S''*(Aa;,u). However, the second 
expression for 5'''(Aco',f) given in Eq. (45) of Ref. [2] is only valid for Vco > Vsi- For 
finite values of U/t this excludes electronic densities in the vicinity of 1. 

Generalization of 5''*(Aa;,f) for the whole range of electronic densities of the 
metallic phase, including densities such that the spin and charge velocities obey the 
inequality Vsi > fco, leads to the following expression, 

Bi,i(^Au;,v) = '-^ I duj' I dk' 

X 



J —sg-a(v)/S.ijj /va 



B^. (^^/^ - k', Aa; - cv') B^^^^ (k\ j) . (23) 
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Figure 4. The value of the cO and si branch hnes exponents of the one-electron 
removal weight distribution as a function of momentum for n = 0.59, m —^ 0, and 
several magnitudes of U/t. The exponents of index c and s refer to the cO and si 
branch lines, respectively. 



In this equation, in those provided in Appendix A, and in the whole our analysis until 
subsection 3.2 we use a notation for the charge branch cO and spin branch si such 
that au = cO, si when the branch index ai? defined by Eq. (Q reads cO = si, cO, 
respectively. Moreover, in equation (123|) / = +1 or / = — 1 according to the A/'-electron 
spectral function ([2]) under consideration, Auj = {uj — lAE), Ak = {k — lAP), and lAE 
and lAP correspond to the general energy and momentum spectra given in Eqs. (28) 
and (29) of Ref . [2] , respectively, which is generated by the elementary processes (A) and 
(B), i = 0, 1, 2, and for i > the index i = 1, 2, ... is a positive integer number which 
increases for increasing values of the number of extra pairs of creation and annihilation 
rotated-electron operators in the expressions of the pseudofermion operators associated 
with the function 5''*(Au;, v) relative to that of the pseudofermion operators of the i = 
function &'^{Auj^v) [2j. Expression ( l23l) is valid for the whole (fc, u;)-plane, except for 
k and uj values such that uj ~ LVav{k — Ik^) + Iuoq, where av = cO, si, l = ±1, ujq and 
fco are provided in Eqs. (32) and (34) of Ref. [2j . 

The velocity v appearing in the argument of the function fl23|) plays an important 
role in our study and is given by, 

Au (u - lAE) , , , ^ , , , 

^ = ^ = 1^—^^^; sgn(w) 1 = sgn(AA;) / ; \v\>Vo,u- (24) 

The inequality |t>| > v^u follows from the structure of the spectrum of Eq. (31) of Ref. 

When at least one of the two parameters 2A^3 ^^^d two parameters 2A^p is finite, 
use of the general expression (123|) for the function &''^{Auo^v) for small finite values of 
lAu = l{uj — lAE) with B^Q^^{k' ,uj') and BQ__{k' ,uj') equaling the suitable expressions 
provided in Eqs. f|T9l) . fl20l) . and fl2T|) leads to the following asymptotic behaviour for 
that convolution function. 



B'^\Auj, v) ^ , ^v^^ 9 /Acu ; . (25) 

Here i = 0,1,2,..., / = ±1, and (o is the functional provided in Eq. (fTOj) and the 
pre-factor function Fq{z) is given in Appendix A. In turn, if all four above parameters 
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vanish we find by use of expression fl23l) with i?g^ (A;', tu') and BQ__{k' ,uj') of the form 
given in Eq. ( l2Ti) that, 

&^'{Auj, v) = B^^\Auj, Ak) = 2nDo 5{Aiu)6{Ak) , (26) 

where Dq is defined in Eq. ( !T4l) . in agreement with the first expression of Eq. (59) of 
Ref. [2] for Co = 0. 

Depending on the values of the two parameters 2A'^^ and two parameters 2Af^, 
we consider in Appendix A seven cases where by use of Eqs. ( |T9l) . ( l20l) . f l2Tl) . ( l23i) . and 
(!25|) a set of alternative expressions for the pre-factor function Fo{z) on the right-hand 
side of the latter equation can be derived and are given in that Appendix . Expressions 
flA.ip - flA.7l) of the same Appendix provide a generalization of the pre-factor function 
Fo{z) given in Eq. (62) of Ref. |2j. The latter function corresponds to the general 
function ( lA.ll) of Appendix A for au = cO [25]. The recent studies of Ref. [1] on the 
finite-energy spectral features of the organic compound TTF-TCNQ use the functions 
(IA.ll) -( )Xr7ll of that Appendix derived in this paper. 

Importantly, general spectral-function expressions for all electronic densities of the 
metallic phase are obtained by replacing the spin velocity Vsi by Vap in the limits of the 
variable z integrations of expressions (66) and (B.14) of Ref. [2] and both in the limits 
of the variable z integrations and arguments of the theta-functions of expressions (71), 
(B.17), and (B.18) of the same reference. After such replacements, Eqs. (66), (68), 
(70), and (71) of Ref. [2j with the pre-factor function Fo{z) given in Eqs. flATTjl - flATTl) 
of Appendix A provide general spectral-function expressions for all electronic densities 
of the metallic phase and spin densities in the range < m < n. However, note that 
while the branch index au is that defined by Eq. (Q, the branch indices au and a'u' 
of Eqs. (66) and (B.14) of Ref. [2] correspond to the two created quantum objects and 
in contrast to the branch index au of Eqs. (I23l) - fl26l) and flA.l|) - flA.7l) of Appendix A 
there are no restrictions imposing that such branch indices are different or equal to the 
branch index au, the same applying to the branch index au of Eqs. (70), (71), (B.17), 
and (B.18) of that reference. 

It is straightforward to show that the double Fourier transform of the function 
ffTSjl equals that of the function f|T71) as 2A^^ ^0. It follows that the asymptotic 
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Figure 6. Full PDT one-electron removal spectral- weight distribution associated 
with the cO and si momentum deviations of Eqs. (|32p and (j33p and corresponding 
functionals given in Eqs. p4|) and ((35| for U/t = 100, n = 0.59, and m ^ 0. 



expression of the correlation function given in Eq. (32) of Ref. [7j is of the general form 
provided in Eq. (52) of that reference independently on whether all four parameters 
2A^^ where au = cO, si and l = ±1 are finite or some of these parameters vanish. 
Thus, the different expressions flA.ll) - flA.7|) of Appendix A for the pre-factor Fq{z) 
of the convolution function ( l25i) are not associated with different pre-factors for the 
corresponding correlation-function expression (52) of Ref. [7J. The pre-factor of the 
latter function is always of the form given in Eq. (50) of that reference |25j . 



3.2. Explicit spectral-function expressions in the vicinity of the border lines 

The singular features of the A/'-electron spectral functions (I2l) are of power law type. 
The power-law branch-line singular features of general form given in Eqs. (68) and (71) 
of Ref. [2] are controlled by non-classic interaction, density, and momentum dependent 
exponents. The spectral feature of general form given in Eqs. (66) of that reference 
can also include power-law singular features called border lines. However, the studies 
of Ref. [2] did not provide an explicit general expression for the k and u dependence in 
the vicinity the border lines. Here we provide such an expression and find that for the 
border lines the power-law exponent has an universal value given by —1/2. 

The shape u = ujBL{k) of a border line is defined by the following parametric 
equations [2], 

k =lka'y',a"u"{q'-, q")^v^,,j{q'),v^„^„{q")- (27) 

Here q' and q" are the bare-momentum values of the a'v' and a'V", respectively. 
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Figure 7. The same PDT one-electron removal spectral-weight distribution as in Fig. 
6 for the values U/t = 4.9, n = 0.59, and m ^ suitable to the TCNQ related spectral 
features. 



pseudofermion or pseudofermion-hole scattering centers created by the processes (A), 
c'l = +1 for pseudofermion creation and c[, c" = — 1 for pseudofermion-hole creation, 
uq is the finite-energy parameter given in Eq. (32) of Ref. [2], and ka'u',a"u"{(l' , <?") is the 
momentum spectrum provided in Eqs. (64) and (65) of that reference. As for Eqs. (66) 
and (B.14) of Ref. [2] and in contrast to the branch index au of Eqs. fl23l) - flA.7l) . there 
are no restrictions imposing that the branch indices a'u' and a'V" of the two created 
scattering centers are different or equal to the branch index ai? defined by Eq. (Q. 

By use of methods similar to those used in Ref. |2] for other weight distributions, 
we find that the spectral function has the following singular behaviour in the vicinity 
and just below (/ = -|-1) or above (/ = —1) the border line. 



2e(n -i[ujBL{k) 



UJ 



TT Cc Cs Co{q', q") \/\aa'u'{,q') \ + \aa"u"{q")\ 



X 



X 



dz Fq{z) 



Co(q',q") f2l[ujBLik) ~uj] 



VcO Vsl 



(28) 



where aav{q) = dvav{q) / dq. In turn, just above (/ = -|-1) or below (/ = —1) the line the 
weight distribution reads. 



(fc, Uj) 



2&(n-l[uj-ujBL{k)] 



7r Cc Cs Co{q', q") \/|aaV'(g')l + (<?")! 
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Figure 8. The region of the (k,uj) plane with a finite one-electron addition spectral 
weight obtained by running over all the momentum values of the cO and si momentum 
distribution deviations of Eqs. ([32| and (|33p for U/t = 100, n = 0.59, and m ^ 0. 
The corresponding branch and border lines are also shown. 



X <^ Q[Vs,p [1 - -] - \Va'u'{q)\ 



n 



1_ 

cti> 



dz Fq[z) 



X 



sgn(<;') 1 



Va'u'{q') \ - Vau [l 

dz Fq{z) 



1[UJ - LOBLjk)] - 



Co(g',g") 



2l[u - UJBL{k)] 



1[UJ - UJBL{k)] 



Co{g',g") 



VcO Vsl 



(29) 



where 6{x) = for x < and 6{x) = 1 for x > 0, i7 is the energy cut off of the 
elementary processes (C) defined in Ref. [2], = Coil': o") is the functional given in 
Eq. ffTOl) for the excited states which contribute to the weight distribution (66) of that 
reference, and Cc and Cg are defined in Eq. (68) of Ref. [3]. 



4. APPLICATIONS OF THE PSEUD OFERMION DYNAMICAL 
THEORY TO THE SPECTRUM OF TTF-TCNQ 

As discussed in Ref. jl], the spectral- weight distribution of TTF-TCNQ is fully 
determined by the occupancy configurations of the c and si pseudofermions in the 
one-electron excited states. The studies of that reference reveal that for electronic 
density n = 1.41 the electron removal spectrum calculated for t = 0.35 eV and 
U = 1.96 eV {U/t = 5.61) yields the best agreement with the TTF related experimental 
dispersions. In turn, for electronic density n = 0.59 an almost perfect agreement with 
the TCNQ related experimental dispersions is reached for the finite-energy-electron- 
removal spectrum calculated for t = 0.40 eV and U = 1.96 eV {U/t = 4.90) [U [18] . 

If one profits from the model particle-hole symmetry, one can achieve the one- 
electron removal spectrum at n = 1.41 and a given U/t value from the one-electron 
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Figure 9. The same region of ttie (fc, us) plane and branch and border lines as in Fig. 
8 for U/t = 5.61, n = 0.59, and to 0. 

addition spectrum at n = 0.59 and the same value of U/t. The full theoretical line 
shape of the one-electron removal spectral-weight distribution found in Ref. [1] to fit 
the corresponding spectral features of TTF-TCNQ is shown in Fig. 1. Note that here 
such a spectrum is rotated relative to that shown in Fig. 2 of that reference. 

Let us consider the processes that contribute to the spectrum of Fig. 1. Each of the 
corresponding spectral-weight contributions to the full one-electron spectrum plotted in 
that figure are obtained by use of the PDT expressions provided in the previous section 
and Ref. [2]. In addition to the values of U/t and n specific to TTF-TCNQ, some of 
the spectral contributions plotted below refer to variations in these parameters. That 
allows us to study how the spectral features evolve under such variations. The study 
of the spectral-weight peaces contributing to the full spectrum plotted in Fig. 1 allows 
the identification of the contributions from the different expressions considered in the 
previous section and Ref. [2] , including the border line expressions derived in this paper 
and provided in Eqs. fl28|) and fl29l) . 

The one-electron removal spectral-weight distribution of Fig. 1 has two main 
contributions referring to the TCNQ and TTF related spectral features, respectively. 
The TCNQ related spectrum is obtained by considering the one-electron removal 
spectral function for n = 0.59; t = 0.40eV;f//t = 4.90. In turn, we profit from 
the model particle-hole symmetry to derive the TTF related photoemission spectrum 
for n = 1.41; t = 0.35 eV;U/t = 5.61 from the corresponding one-electron addition 
spectrum for n = 0.59; t = 0.35 eV; U/t = 5.61 and changing the energy and momentum 
signs. 

Nearly the whole one-electron removal spectral weight corresponds to excitations 
described by the following deviations from the ground-state N^i, cxv pseudofermion 
numbers and N^^ av pseudofermion-hole numbers where av = cO, si, 

AiVeo = -AN':, = -1 ; A^^si = -AN^i = -1 , (30) 

and with the numbers Qauf^ appearing on the right-hand side of Eq. ([7]) reading, 

Q%/2 = ±n/2; Q'J2 = . (31) 
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Figure 10. The value of the cO and si branch lines exponents of the one-electron 
addition spectral-weight distribution as a function of momentum for n = 0.59, m —^ 0, 
and several magnitudes of U jt. The exponents of index c and s refer to the cO and si 
branch lines, respectively. 



If we consider bare-momentum continuum values, the general deviations read, 

Lj 

27r 

AiV,i(g) = - — 5(g - g',) ; g\ G [-A;^^, +fc^J , (32) 

and thus 

27r 

AiV,i(g) = -— 5(g - g\) ; g\ G \-kp, (33) 

as m — i> for the initial ground state. Such deviations and numbers (13T!) are then 
used in the general functional Qav{(lj)l'^ given in Eq. (jTj), whose functional Q'^^{qj)/2 
is defined in Eq. (jl]). Such a procedure leads to, 

Qco(g)/2 = ±7r/2-7r<l>co,co(g,gi) - 7r<l>co,co(g,T2A;F)/2 

+ 7r$co,co(g, ±2kF)/2 - 7r$eo,si(g, g'l) , (34) 

and 

Q,i(g)/2 = - 7r<l>,i,co(g, gi) - vr$si,co(g, T2A;f)/2 

+ 7r$,i,co(g, ±2A;i.)/2 - 7r<l>,i,,i(g, q\) , (35) 

respectively. 

The creation of the holes in the cO and si bands whose energy dispersions are 
plotted in Figs. 6 and 7 of Ref. [15], respectively, corresponds for U /t = 100, n = 0.59, 
and m — > to a finite spectral-weight distribution in the region of the {k, u) plane 
shown in Fig. 2. The lines associated with the branch lines and border lines in the 
finite-weight distribution are also shown. Most lines of the one-electron spectral- weight 
distribution plotted in Fig. 1 corresponding to power-law singularities are of the branch 
line type. In contrast to the border-line singularities of expressions (1281) and (|29i) . whose 
exponent —1/2 is momentum and U/t independent, that controlling the branch-line 
power-law singularities is both momentum and interaction dependent. Near a branch 
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Figure 11. The same exponents as in Fig. 10 as a function of momentum for U /t = 10, 
TO ^ 0, and several values of n. 



line the general expression of the weight distribution is that provided in Eq. (70) of Ref. 
[2j. Only for momentum and interaction values where the exponent of that expression is 
negative is the corresponding line called a branch line. For the limit m — >^ considered 
here the pre-factor function -^0(2;) appearing in such a branch-line expression is that given 
in Eqs. flA.ll) - flA.7p of Appendix A. The expressions provided in the latter equations 
are a generalization of the pre-factor function Fq{z) of Eq. (62) of Ref. [2]. The latter 
applies solely when < n < 1 and < m < n. 

In figure 3 the same finite one-electron removal spectral-weight region and 
corresponding lines are shown for the U/t = 4.9, n = 0.59, and m values suitable 
to the TCNQ related spectral features. Consistently with the U/t dependence of the 
energy bandwidth of the si pseudofermion dispersion plotted in Fig. 7 of Ref. [15], note 
that the border line connecting the minimum energy points of the two cO-branch lines 
has Vco{q) = Vsi{q') ~ 0. Upon decreasing the U/t value to that U/t = 4.9 suitable to 
the TCNQ related spectral features of Fig. 1, such a border line acquires a curvature. 
The curvature achieved at U/t = 4.9 is that behind the agreement found in Ref. [1] 
between the PDT predictions and the photoemission spectral features of TTF-TCNQ. 
That agreement is not obtained for larger or smaller values of U/t. 

The (negative-exponent) cO branch line of the weight distributions of Figs. 2 and 

3 runs between the excitation momenta —kp and Skp and the si branch line from —kp 
to kp. The corresponding power-law momentum dependent exponents are plotted in 
Fig. 4 for n = 0.59, m — 0, and several values of U/t. The cO branch line segment 
between —kp and is folded in the positive momentum region. While the negative cO 
branch line exponent is smaller for U/t large, the negative si branch line exponent is 
smaller for smaller values of U/t. In order to illustrate how systematic variations in the 
parameters lead to the evolution of the spectral features, the same exponents as in Fig. 

4 are plotted in Fig. 5 as a function of momentum for f//t = 10, m ^ 0, and several 
values of n. 

The use in the general PDT expressions provided in this paper and in Ref. ^ of 
the specific cO and si momentum deviations of Eqs. fl32|) and fl33l) and corresponding 
functionals given in Eqs. (!34l) and ( !35l) leads for U/t = 100, n = 0.59, and m — to 
the one-electron removal spectral- weight distribution plotted in Fig. 6. Such a large- 
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Figure 12. Full PDT one-electron addition spectral-weight distribution associated 
with the cO and si momentum deviations of Eqs. (I38p and ()39|1 and corresponding 
functional given in Eqs. (gO]) and (gl]) for U/t = 100, n = 0.59, and m ^ 0. 

U/t distribution is quite similar to the function B{k,u!) plotted in Fig. 1 of Ref. [I6] 
for U/t ^ oo. However, the method used in the studies of the latter reference does 
not apply to finite values of U/t. The one-electron removal spectral- weight distribution 
related to the TCNQ spectral features is plotted in Fig. 7 and refers instead to f//t = 4.9, 
n = 0.59, and m 0. It is a part of the full one-electron spectral weight distribution 
plotted in Fig. 1. 

Most singular behaviours of the one-electron removal spectral-weight distributions 
of Figs. 6 and 7 refer to cO and si branch lines whose negative power-law exponents 
magnitude depends on U/t and the momentum. These exponents are plotted in Figs. 
4 and 5. In turn, the weaker lines connecting in Figs. 2 and 3 the minimum energy 
points of the two cO lines are border lines. The singularities of the latter lines correspond 
instead to the U/t and momentum independent exponent —1/2, as given in the general 
border-line expressions (l28l) and (l29l) derived in this paper. 

Next let us consider the lower-Hubbard-band one-electron addition spectrum 
for n = 0.59; t = 0.35 eV; U/t = 5.61 from which one derives the TTF related 
photoemission one-electron removal spectrum for n = 1.41; t = 0.35 eV;U/t = 5.61. 
Nearly the whole lower-Hubbard-band one-electron addition spectral weight corresponds 
to the excitations leading to the following deviations from the ground-state Nau 
pseudofermion numbers and Nj^^ av pseudofermion-hole numbers where av = cO, si, 

AN,o = -AN^Q = +l; ANsi = 0; AN^, = +1 , (36) 
and with the numbers Q^^tv/I appearing on the right-hand side of Eq. (I7j) reading, 

Qy2 = 0; Q°i/2 = ±7r/2. (37) 
If we choose the excitation branches that survive for m ^ 0, the general deviations are 
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Figure 13. The same PDT one-electron addition spectral-weight distribution as in 
Fig. 8 for the values U/t = 5.61, n = 0.59, and m suitable to the TTF related 
spectral features. 



for bare-momentum continuum values given by, 

AAr,o(g) = —6{q- qi) ; qi G [-2kF, +n] and qi G [-vr, +2^^] 
AiV,i(g) = - 5(g - q\) + j5iq + kp^) 

+ ^ 5(g - kpi) ; q\ G [-kp^, +kp^] , (38) 

and thus 

AiVsi(g) = - ^5{q- q\) + ^5{q- kp) 

+ ^5{q + kp) ■ q\ G [-kp, +kp\ (39) 

as m — > for the initial ground state. Alike for the above one-electron removal case, such 
deviations and the numbers fl371) are used in the general functional Qauilj)/'^ provided 
in Eq. ([7]) and its functional QtuiQj)/'^ defined in Eq. (j4]). One then finds, 

<5co(g)/2 = n^co,co{q,qi) - vr<l>co,si(g, g'l) 

+ vr<l>eo,.i(g, -kFi)l2 + 7r$,o,si(g, +/i;Fi)/2 , (40) 

and 

Q,i(g)/2 = ± 7r/2 + 7r$,i,co(g, gi) - 7r$,i,,i(g, g'^) 

+ vr<l>,i,,i(g, -A;Fi)/2 + 7r$,i,,i(g, +A;Fi)/2 , (41) 

respectively. 

The creation of the cO pseudofermion and si pseudofermion hole refers for f//t = 
100, n = 0.59, and m — > to a finite spectral-weight distribution in the region of 
the {k,uj) plane shown in Fig. 8. The corresponding lines associated with potential 
branch lines and border lines are also plotted. In Fig. 9 the same finite spectral- weight 
region is shown for the U/t = 5.61, n = 0.59, and m — values suitable to the TTF 
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related spectral features. The border line connecting the maximum energy point of 
the cO-branch line at momentum 2kF with the maximum energy point of the finite 
spectral- weight distribution at momentum vr has in Fig. 8 Vco{q) = Vsi{q') ~ 0. Upon 
decreasing U/t to the value U/t = 5.61 suitable to the TTF related spectral features of 
Fig. 1 that border line acquires a curvature, as shown in Fig. 9. Such a curvature and 
corresponding value of U/ 1 are those behind the agreement achieved in Ref . ^ between 
the PDT predictions and the photoemission spectral features of TTF-TCNQ. 

The cO branch lines of the weight distributions of Figs. 8 and 9 run between 
the excitation momenta kp and vr and between Skp and vr, respectively. In turn, the 
maximum extension of the si branch line is achieved for large values of U/t when it runs 
from kp to near Skp. The corresponding power-law momentum dependent exponents 
are plotted in Fig. 10 for n = 0.59, m — 0, and several values of U/t. While the cO 
branch lines exponents are negative, the si branch line exponent is negative for momenta 
smaller than approximately kp and 2kp for U/t = 5.61 and U/t = 100, respectively. 
Again in order to illustrate the evolution of the spectral features upon varying the 
parameters, the same exponents as in Fig. 10 are plotted in Fig. 11 as a function of 
momentum for U/t = 10, m ^ 0, and several values of n. 

The general PDT spectral-function expressions provided in this paper and in Ref. 
[2] lead for U/t = 100, n = 0.59, and m — ^ and the specific cO and si momentum 
deviations of Eqs. (!38l) and ( |39l) and corresponding functionals provided in Eqs. (HOll 
and (HTj) to the one-electron addition spectral-weight distribution plotted in Fig. 12. 
The large-f//t distribution plotted in that figure is quite similar to the function A{k,u!) 
plotted in Fig. 1 of Ref. [16] for U/t — > oo. The one-electron addition spectral- 
weight distribution related to the TTF spectral features is plotted in Fig. 13 and refers 
instead to U/t = 5.61, n = 0.59, and m —>■ 0. After a straightforward particle-hole 
transformation, the latter distribution leads to the one-electron removal spectral-weight 
distribution suitable to the TTF related spectral features. The latter corresponds to 
U/t = 5.61, n = 1.41, and m — and is part of the full one-electron removal spectral- 
weight distribution plotted in Fig. 1. 

Most singular behaviours of the one-electron addition spectral-weight distributions 
of Figs. 12 and 13 correspond to cO and si branch lines whose negative power-law 
exponents are dependent on the U/t and momentum values. These exponents are plotted 
in Figs. 10 and 11 and the power-law singularities occur for the momenta where they 
are negative. In turn, the line connecting in Figs. 8 and 9 the maximum energy point 
of the cO-branch line at momentum 2k p with the maximum energy point of the finite 
spectral-weight distribution at momentum vr is a border line. The singularities at the 
latter line correspond to the U/t and momentum independent exponent —1/2 and near 
it the line shape is of the general form provided in expressions (l28l) and (!29l) . 

The four functionals 2A^^ where au = cO, si and t = ±1 given in Eq. (|TOl) are in the 
case of the one-electron excitations behind the spectral-weight distribution plotted in 
Fig. 1 fully controlled by the functionals provided in Eqs. ( IMI) and ( l35i) for U/t = 4.90, 
n = 0.59, and m ^ for the TCNQ related spectral features and in Eqs. (j40l) and 
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(HTj) for U/t = 5.61, n = 0.59, and m — for the TTF related spectral features. Since 
the magnitude of the latter functionals is a function of two momenta qi and q\ of the 
deviations of Eqs. (l32l) and ( l33l) or Eqs. (138|) and ( 139|) belonging to the cO and si band, 
respectively, it is different for each point of the {k, u) plane. Moreover, for some of the 
latter points there are contributions from more than one pair of momenta of qi and q\. 

We emphasize that depending on whether none, one, or two of the two above 
momenta qi and q'l is or are Fermi points, the corresponding functionals 2AJ;^^ have 
different form and the excitation contributes to line shapes described by the expressions 
provided in Eqs. (66), (70), and (68) of Ref. [2], respectively. All such expressions 
involve the pre-factor function Fq{z) given in Eqs. (lA.ll) -( IX]7l) of Appendix A. The 
expressions provided in these equations are a generalization of the pre-factor function 
Fo{z) provided in Eq. (62) of Ref. [2]. The latter applies when the four functionals 
2A^^ where au = cO, si and l = ±1 are finite. However, in the limit m — that the 
spectral-weight distributions studied in this paper refer to one must use in the general 
expressions provided in Eqs. (66), (70), and (68) of Ref. [2] the pre-factor function 
Fo{z) of Eqs. flCT- flAJ]) of Appendix A. 

Note that the above mentioned general expressions provided in Eqs. (68) and (70) 
of Ref. [2] describe the line shape near well-define points and lines in the {k, u) plane, 
respectively. When the corresponding exponents are negative, such spectral features 
refer to the point and branch-line singularities of the weight distributions plotted in 
Figs. 1, 6, 7, 12, and 13. In turn, the only type of singularity occurring within the 
line shape described by the above mentioned general expression provided in Eq. (66) 
of Ref. [2] is that associated with the border lines studied in the previous section. The 
corresponding expressions ( !28l) and (!29ll refer to the vicinity of such border lines. 

The complementary studies of Ref. [18] provide the specific expressions of the four 
functionals 2A^^^ where au = cO, si and t = ±1 given in Eq. ffTOl) for each of the 
point and branch-line singularities of the TCNQ related one-electron removal spectral 
features. The functionals given in Eqs. (HO!) and (HTl) are the basis of the derivation of 
similar expressions for the TTF related spectral features. Those are straightforwardly 
derived by choosing particular pairs of momentum values qi and q'l in the cO and si 
momentum deviations of Eqs. ( !38|) and ( !39l) . respectively. As discussed above, when 
such choices involve none, one, or both such momentum pairs being a cO or si Fermi 
point, the corresponding four functionals 2A''^^ have different forms and the excitation 
contributes to hne shapes described by the expressions provided in Eqs. (66), (70), and 
(68) of Ref. [2] for the two-dimensional weight distribution, vicinity of branch lines, and 
vicinity of singular points, respectively. 

Finally, in Fig. 14 we plot the theoretical lines corresponding to the sharpest 
spectral features considered in Fig. 1 but omit the corresponding detailed spectral- 
weight distribution over the {k, a;)-plane provided in that figure. The figure also 
displays the experimental dispersions in the electron removal spectrum of TTF-TCNQ 
as measured by ARPES in Ref. [9]. The border line connecting the maximum energy 
point of the cO-branch line at momentum 2k p with the maximum energy point of the 
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Figure 14. Experimental peak dispersions (grey scale) obtained by ARPES on TTF- 
TCNQ along the easy-transport axis as given in Fig. 7 of Ref. [S] and matching 
theoretical branch and border lines. (The Z-point corresponds to the momentum 
k = TT.) The corresponding detailed theoretical spectral-weight distributions over the 
whole {k, tt')-plane are plotted above in Fig. 1. While the theoretical charge-c" and 
spin-s" branch lines and c — s border line refer upon the particle- hole transformation to 
the TTF related spectral features of Figs. 9 and 13, the charge-c, spin-s, and charge-c' 
branch lines correspond to the TCNQ related dispersions of Figs. 3 and 7. 



finite spectral-weight distribution at momentum vr in Fig. 9 leads upon the particle-hole 
transformation to that called c — s line in Fig. 14. In turn the border line connecting 
the minimum energy points of the two cO lines of Fig. 3 associated with the TCNQ 
related spectral features is not marked in Fig. 14, yet the corresponding experimental 
spectral weight is clearly visible. 

5. DISCUSSION AND CONCLUDING REMARKS 

In this paper we have generalized the closed-form analytical expressions for the finite- 
energy one- and two-electron spectral-weight distributions of a ID correlated metal 
with on-site electronic repulsion introduced in Ref. [2] to all electronic densities of 
the metallic phase and zero spin density. Moreover, we have studied the particular 
form of the expressions derived here for the processes contributing to the one-electron 
spectral-weight distributions related to the TTF and TCNQ stacks of molecules, 
respectively, in the quasi-lD organic compound TTF-TCNQ investigated in Ref. [1]. 
The corresponding full theoretical one-electron spectral-weight distribution plotted in 
Fig. 1 agrees quantitatively for the whole experimental energy bandwidth with the 
observed one-electron spectral features shown in Fig. 14. 

Other applications of our finite-energy spectral-weight-distribution expressions to 
several materials, correlated quantum systems, and spectral functions are in progress. 
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This includes the use of our theoretical results in the two-atom spectral-weight 
distributions measured in ID optical lattices. Such studies also involve the Hubbard 
model but with the electrons replaced by fermionic spin- 1/2 atoms. While for the 
one-electron features the branch lines play the major role, in the case of two-electron 
spectral functions such as the charge dynamical structure factor the two-pseudofermion 
spectral features of the form given in Eq. (66) of Ref. |2j lead to the main contributions. 
Generalization of that equation to all electronic densities of the metallic phase and to 
zero spin density involves the use of the following expression, 



dz Fq{z) 



{ " y 



■ l = ±l, (42) 



Co(g', q") \Va'u'{q') - Va"u"{q")\ 

with the pre-factor function Fq{z) given by the expressions (lA.ip - (lA.7p of Appendix A 
introduced in this paper. Alike in Eq. fl28|) . here q' and q" stand for the bare-momentum 
values of the created aV and a'V" pseudofermions or holes, respectively. Indeed, the 
singular border-line function (!28l) is a particular case of the general function provided in 
Eq. (1421) ■ which can be obtained by considering that Va'u'W) = Va"u"{(l") in the latter 
equation. 

For finite values of U/t the dominant term of the two-electron spectral functions 
is often of the general form given in Eq. ( H2l) . However, it can also occur that 
for intermediate values of U/t such functions are the sum of two or three dominant 
functions of the general form (H2l) . The main contributions correspond to the two 
created objects being such that (i) aV = cO and a'V" = si, (ii) aV = a'V" = cO, and 
(iii) a'u' = a"v" = si. The relative importance of the functions (i)-(iii) and whether 
the two created objects are pseudofermions and/or pseudofermion holes depends on 
the specific two-electron spectral function under consideration. Our general expressions 
(IA.ll) -( )Xr7j) of Appendix A for the pre-factor function Ff){z) appearing on the right-hand 
side of Eq. f l42l) allows the evaluation of all two-electron spectral-weight distributions 
for initial ground states with zero spin density, which in many situations is the case of 
physical interest. The branch lines and other spectral features also contribute to the 
two-electron spectral-weight distributions, yet the dominant contributions are of the 
general form given in Eq. (I42l) . 

For instance, the preliminary results of Ref. [20] consider both one- and two-electron 
spectral features, profits from the PDT for the one-chain problem, are consistent with 
the phase diagram observed in the (TMTTF)2X and (TMTSF)2X series of organic 
compounds, and explain the absence of superconductive phases in TTF-TCNQ. The 
studies of that reference combine the PDT for several one- and two-electron spectral 
functions with a Renormalization Group analysis to study the instabilities of a system 
of weakly coupled Hubbard chains. For low values of the onsite repulsion U and of the 
doping 6 = (1 — n), the leading instability is towards a superconducting state. The 
process includes excited states above a small correlation pseudogap. Similar features 
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Figure 15. Temperature-pressure phase diagram of the (TMTTF)2X and 
(TMTSF)2X series of compounds. It includes non- Fermi hquid (NFL), insulating 
Mott-Hubbard (Mott-CO), spin-density- wave (SDW), spin-Peierls, Fermi-liquid (FL), 
and superconducting (SC) phases. 



appear in extended Hubbard models in the vicinity of commensurate fillings. The 
theoretical predictions of such studies are consistent with the phase diagram observed 
in the (TMTTF)2X and (TMTSF)2X series of organic compounds represented in Fig. 
15. From the results of Ref. [20] one can infer that the low-temperature phase of 
the coupled-chain system will show long-range superconducting order. However, the 
precise nature of this phase, and the symmetry of the order parameter is dependent 
on the arrangement of the chains within the material. The investigations of Ref. [20] 
use the exponents obtained from the PDT for the one-chain problem and confirm that 
in ID non-Fermi liquids weak inter-chain hopping can induce superconductivity. At 
low temperature these materials show a spin-Peierls or spin-density-wave phase. Under 
pressure, the (TMTSF)2X compounds are driven to a superconducting phase, which is 
removed again if one further increases the pressure. 

In addition to the finite-energy spectral features of the organic compound TTF- 
TCNQ [1] and the preliminary studies of Ref. [20] on the (TMTTF)2X and (TMTSF)2X 
series of organic compounds, the general spectral-function expressions derived in this 
paper will be used elsewhere in the study of specific one- and two-electron spectral 
functions and its quantitative application to the unusual spectral properties of other 
low-dimensional materials and systems. While the studies of this paper considered the 
ID Hubbard model, which describes successfully some of the exotic properties observed 
in low-dimensional materials [9], [18], [19], [20] , our results are of general nature for many 
integrable interacting problems [1] and therefore have wide applicability. 
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Appendix A. The pre-factor function Fq{z) 

In this Appendix a set of alternative expressions for the pre-factor function Fq{z) on 
the right-hand side of Eq. ([251) derived by use of Eqs. (EQl), (EI]), (|23l). and (|25l) 

according to the values of the two parameters 2A^3 and two parameters 2A^p are given. 

(i) When the four parameters 2A^3 where au = cO, si and l = ±1 are finite we find, 



Fo{z) = 2DoJ— / dx 

V Vau Jo 

6 ( 1 — a; + sgn(z) l' 



dy 



n 

l'=±1 

9 + sgn(z) l' y 



Vapz -f^y 



1 — a; -|- sgnf^;) l' 



Vn 



X + sgn(z) L y 



2A'„-1 



(ii) When 2A^p = and the remaining three parameters are finite. 



Fo{z) = 2Do^ 



dy 



'Vay J -I 

e ( I [sgn(2)?/ -Vo^^iz- [zvc,u - sgT).{z)y] 



Vn 



21 



z — sgnU 



r(2A 

V 



n 

t'=±i 



Vqiv 

efl + sgn(z) 



y - LVauZ 



r(2A^j 



(A.l) 



1 + sgn(2;) 



/ , — Vqii/ 

i + L 



y- iv^pz 



(iii) When 2A^^ = and the remaining three parameters are finite. 



Fn(z)= 2Dn 



-1 'diLSgn{z)y 
dy ^ 



(A.2) 



n 

r'=±i 



6(1 — sgn(z) 



r(2A^, 



— LSgn{z)2y 
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r(2AL 



y + l'vau z^ 
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X 



(iv) When both 2A^p — and the remaining two parameters are finite, 



1 — sgn(2;) 



L + V 



y + LVapZ 



(A.3) 



Fo{z) = Do.I^q(—-\z\) 

V ^Vau f 



X 



n 



'OLV 



1 + l!Vo,v Z 



^Al^r(2A^j vv 

(v) When both 2A^^ = and the remaining two parameters are finite. 



(A.4) 



Fo(^)=Do./— n 



-tz 



1- LVn,iyZ 



(vi) When 2A~^ = 2A^^ = and the remaining two parameters are finite, 

2Do I v„„ - LiVi 

F^Kz) = 



0.(A.5) 



X 



1-Co 



r(2A^,) r(2A^ 



av I 
2A^--1 



iZ\ 



(A.6) 



(vii) When both 2A^^ > (and 2A^p > 0) and the remaining three parameters vanish, 

2A'„-1 



S z 



(A.7) 



V ^au J \ 

(and a similar expression with av^i^ and av replaced by aP, i, and av^ respectively). 
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